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I. INTRODUCTION 

Superconducting circuits are a promising approach 
to building a large-scale quantum information proces- 
sor. Over the past ten years, quantum coherence times 
have improved by two orders magnitude, from ns to fis 
timescales [l|, H, y, 13, H, ■ With this improvement has 
come increased attention to the fundamental quantum 
processes that arise when these circuits are controlled 
by microwave fields. A recurring theme in recent ex- 
periments has been to characterize the multiple quan- 
tum levels that can be excited, in either the frequency 
or time domain. On the one hand, these extra levels can 
interfere with ideal qubit operation. There have been 
many theoretical studies of the imperfections that arise 
due to higher energy levels, a phenomenon called "leak- 
age" On the other hand these higher levels can also 
be used advantageously, either to mediate quantum inter- 
actions between qubits \§\ or to process quantum infor- 
mation with higher dimensional quantum systems called 
qudits . Recently, multiple levels of a superconducting 
phase qubit have been addressed by multi-frequency con- 
trol fields to emulate a quantum spin (with spin > 1/2) 
(lo| . From either perspective, it is an important task to 
develop theoretical tools to model these quantum pro- 
cesses simply and accurately. 

Most theoretical work has focused on the deviations 
from ideal qubit behavior during Rabi oscillations, and 
how these can be mitigated by pulse shaping techniques 
[Til . [T2[ . Optimal control theory has also been applied 
to this problem [l^ [3, [l^, and recent work has indi- 
cated that arbitrarily fast control is possible using cer- 
tain choices of pulses [l^ . The presence of higher levels 
is also problematic for coupled qubit operation. These 
arose in the study of coupled phase qubits: the spec- 
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troscopic signatures were analyzed in [17|, while a non- 
adiabatic controlled phase gate using the higher levels 
was first proposed in [1]. 

Experimentally, the effect of the higher levels in a su- 
perconducting circuit has been demonstrated in trans- 
mon circuits .18], both in single-qubit operations 19] and 
recently in a two-qubit controUed-phase gate [20j similar 
the phase qubit gate described above. For phase qubits, 
multilevel Rabi oscillations [111, and multi-photon 
Rabi oscillations 23] have been analyzed in some detail, 
while a sensitive characterization of leakage was demon- 
strated by a Ramsey filter method [IJ. Recently, inter- 
ference effects due to multiple frequency controls were 
demonstrated [1^, realizing effects related analogous to 
electromagnetically-induced transparency [1^, [l^l • 

In this paper we develop a simple theoretical frame- 
work to describe the control of multiple levels in a su- 
perconducting phase qubit using multi-frequency control 
fields. We start from an early proposal to reduce leakage 
during qubit manipulation by resonantly cancelling off- 
resonant transitions to the higher energy levels [2^ . Nu- 
merical simulations are used to demonstrate that this ap- 
proach can optimize a quantum transition on the multi- 
level qubit. These results are explained using the many- 
mode generalization [HI of the Floquet formalism [30| 
for a Hamiltonian that is periodic in time. We show 
that multi-frequency control fields can produce a unique 
quantum interference to optimize the desired transition, 
without complex pulse shaping. We further show how the 
Floquet formalism can describe other interference effects 
when driving multiple transitions. 

This paper is organized as follows. In Section II we 
describe the basic model of a phase qubit. In Section 
III, we introduce the Floquet formalism for a single fre- 
quency control pulse, reproducing the effects that occur 
in three- level Rabi oscillations. This formalism gener- 
alizes the rotating wave approximation, taking a time- 
dependent problem to a time-independent problem (with 
a much larger state space). In Section IV we extend the 
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Floquet formalism to include control fields with multi- 
ple frequencies. This analytical approach is used to opti- 
mize a transition between the first two levels of the phase 
qubit. These ideas are confirmed in Section V through 
numerical optimizations of square and Gaussian control 
pulses. We return to the Floquet formalism in Section VI 
to predict beating effects relevant to the recent spin emu- 
lation experiment [10]. Finally, we conclude our study in 
Section VII, while certain theoretical results are detailed 
in the Appendix. 



The resulting energies and matrix elements, calculated 
using perturbation theory, are found in the Appendix. 
Finally, the driving field has the explicit form 

III. SINGLE-MODE FLOQUET THEORY: 
THREE-LEVEL RABI OSCILLATIONS 



II. PHASE QUBIT HAMILTONIAN 

The phase qubit is generally based on a variation of the 
current-biased Josephson junction [ij. This is described 
by the following Hamiltonian 



H 



Ej{cos-f+{I/Ich). 



(1) 



The dynamical variables are 7, the gauge-invariant phase 
difference, and p-y, its conjugate momentum, subject to 
the commutation relation [7,^7] — ih. The other pa- 
rameters are the junction's bias and critical currents 
/ = I lie and /c, the capacitance C, and the energy scales 
Ej = hlc/2e and E^ = e'^/2C. 

To describe Rabi oscillations, we will let the bias cur- 
rent be time-dependent, of the form / = Idc — Iac{t), and 
restrict the Hamiltonian to the lowest four energy levels 
to find 



H = Hn 



.f{t)X, 



(2) 



where we have divided the Hamiltonian into its unper- 
turbed, time-independent form 



Ho 



£;o 





-Ba 



(3) 



and a set of dimcnsionless matrix elements 



X = 



/ 2^00 2^01 2^02 

Xoi Xll Xl2 Xl3 

X02 X12 X22 X23 

\ Xo3 Xi3 X23 X33 



(4) 



The energy levels En and the matrix elements Xnm can 
be calculated by either diagonalizing the washboard po- 
tential directly, or by some approximation scheme. The 
latter can be efficiently performed by first approximating 
the washboard potential as a cubic oscillator, of the form 



H — tkoQ 



1 2 



1 



Ax' 



whe re hu jQ = ^SE^Ej (l - {Idc/Icf 
l/y/MWs with Ns given by 

AU _ 23/4 
hujo 3 



(5) 



and A = 



N, = 



El 

Er_ 



1/2 



^ Idc 



5/4 



(6) 



For Rabi oscillations in the presence of strong driv- 
ing, there are deviations from two-level behavior that can 
be anal yze d using a three-level model. Previous studies 
[Til [l^. I31I [3^ , using the rotating wave approximation, 
have identified three main features. First, the coherent 
oscillations between the ground and first excited state are 
accompanied by oscillations to the second excited state. 
Second, there is a reduction in a Rabi frequency. Finally, 
there is a Stark shift of the optimal resonance condi- 
tion. All of these effects have been seen experimentally 
[2^ . [2^ . [3^ . In this section we theoretically derive these 
effects by introducting the Floquet formalism [30] . 

First, we let the driving field be given by /(i) = 
Acosujt. Then, we expand the wavefunction as a Fourier 
series 



(8) 



Finally, substituting this scries into the Schrodinger 
equation ihd\9)/dt = with H = Hq + AXcosLut, 

we match terms proportional to e™'^* on each side. The 
resulting equations to be solved are 



ih 



d\ip„ 
dt 



(i/0+riMIV'n> + 7T^^(l^n-l> + IV^n+l». (9) 



Letting |V'n(i)> = e-'-^*/'^|?/'„(0)), we find that these 
coupled equations are equivalent to a time-independent 
Schrodinger equation 7ii?|^') = E\'^) for the infinite state 
1*^) ~ Y^'n=-oo IV'n) ® W) with the Floquet Hamiltonian 
matrix 



F )n.m 



= {Hq + nhuj)5n,'. 



-AX{5 

71 , T71 — 1 



^n.rn+l ) ■ 

(10) 

The labels n and m can be interpreted as photon numbers 
for the driving field, and the overall state as that of the 
combined system and field. 

In general, this approach has replaced a finite- 
dimensional time-dependent problem with an infinite- 
dimensional time-independent problem. To solve the lat- 
ter, we can approximate the infinite matrix by one of 
its sub-blocks. For the problem at hand, the lowest- 
order approximation is to include only three states: |0, 0), 
|1, —1), and |2, —2), where we are using the notation of 
the form |s, n) to indicate the system in state s = 0, 1, 2 
with n = 0, —1, and —2 photons, respectively. Negative 
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photon numbers are allowed here, as these are differences 
from the average photon number in a semi-classical state 
[30l |. After removing an overall constant energy Eq, the 
resulting Floquet matrix takes the form 



floi/2 

f^i2/2 UJ02-'2UJ 



(11) 



where ftfioi = Axqi, ftwoi = Ei—Eq, and hLL!Q2 = E2—Eq. 
For convenience, we also define hu!i2 = E2 — E1; note that 
UJ02 = UJ12+UJ01. Note also that this approach reproduces 
the rotating wave approximation exactly, while includ- 
ing more states allows for systematic corrections due to 
strong rnultiphoton processes, such as the Bloch-Siegert 
shift [30|. The resulting dynamics can be found by di- 
agonalizing the Floquet matrix. For this Hamiltonian 
exact results are available [13, [s^ ; here we will adopt a 
perturbative approach. 

To simplify the following, wc consider the case of near 
resonance with 6 = lu — ujqi <C ^01, and approximate 
fli2 ~ \/2Qoi. For weak driving, the largest scale in the 
problem is A = 2w— u;o2 ~ wqi— W12 ~ 5ll>o/ {36Ns). This 
corresponds to the anharmonicity of the system, being 
inversely proportional to Ng. For Rabi frequencies near 
this value, three-level effects become important 12, 32]. 
Therefore, to see deviations from two-level behavior, we 
use perturbation theory in the small parameters S/flgi 
and fioi/A, starting from the zeroth order eigenstates 
(|0,0) ± |l,-l))/\/2 and |2,-2). Using standard meth- 
ods of perturbation theory, we cornpute the (normalized) 
eigenstates 1^;^) and eigenvalues E^, from which we cal- 
culate the time-dependent amplitudes 



\v,){vi\i;m, (12) 



where we assume that V'(O)) = |OjO). This calculation 
is best done using a computer, as the required order of 
perturbation theory is second order for the wavefunction 
and fourth order for the energy. Alternatively, one can 
expand the exact eigenvalues using the roots of a cubic 
polynomial. In either case, we find that the amplitudes 
satisfy 

ao{t) « cos(m/2) - i sin(m/2) - , (13) 



There are many things to note about this solution. 
First, we observe that both ai(t) and a2{t) are propor- 
tional to sin(ili/2). Thus, transitions from the ground 
to the first excited state leak out to the second excited 
state, with probability 



P2it)^\a2{t)\' 



2A2 



sin2(m/2). 



(17) 



Avoiding this leakage through pulse shapes has been the 
subject of much investigation [Tl], [l^, [3, HBl • In addition 
to this error, however, is the reduction oi pi{t) = \ai{t)\^ 
by the factor depending on (5/rioi- This is due to the 
fact that, when coupled to the second excited state, the 
— > 1 transition is no longer located at S = uj — wqi = 0, 
but rather at S = flQi/{2A), i.e. 



"01 



2(cjoi ~uji2)' 



(18) 



This is the effective ac Stark shift measured in exper- 
iments (23I , |32| . As shown below, it must be compen- 
sated for high-fidelity qubit rotations. Finally, the on- 
resonance Rabi frequency is given by 



~ ^01 1 



"01 



4(cjoi - t^i2)^ 



(19) 



Its reduction is due to the dressed eigenstates of the 
system, and has also been measured experimentally 
21i,i2i[2l,[3i 



IV. TWO-MODE FLOQUET THEORY: 
OPTIMIZED RABI OSCILLATION 

In the presence of a control field of the form 

f{t)^AiCOSUJit + A2COs{LLl2t + <l>), (20) 

the Floquet method can be generalized [2§| to include two 
sets of photon states for the two oscillatory components 
of the field. That is, by performing the double Fourier 
expansion 

00 

!*(<))= E l^m,«.(t))e"^"^*e»^'^^*, (21) 

{ni,n2} = — 00 



ai{t) « -ism{nt/2) 



and 



02 (t) 



02 
"01 

4A2 



6 



-i—^sm{nt/2), 



where the Rabi frequency Q is given by 



~2 



6_ 



^01 
2A 



2A 



(14) 
(15) 

^(h--^) • (16) 



the Schrodinger equation leads to the set of coupled equa- 
tions 

^^ d\i'nin^) ^ [Ho+nihu)i+n2hjJ2)\lpm.n2) (22) 

at 

+ 2^1^(IV'ni-l,n2) + IV'rii-Hl.na)) 

+ iA2X(e^*|V„,,„,-i) + e-^^|V„i,„.+i)). 

This is equivalent to a time-independent Schrodinger 
equation for the infinite state |^) — ns IV'ni,n2) ® 
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where in the following we will assume that |'0(O)) = 
|0,0, 0). The structure of these equations is well de- 
scribed elsewhere [2^ . Here we make the following obser- 
vations. First, to obtain accurate numerical results, one 
must include several photon states in the sum — including 
too few results in a loss of both accuracy and unitarity. 
Second, one can still use perturbation theory to obtain 
useful analytical results, provided one identifies the ap- 
propriate states of the combined system. 
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FIG. 1: Three-level Rabi Oscillation. The probability pi (t) = 
|ai(t)|^ to be in state 1 is shown as a function of time. The 
solid curve is a numerical simulation using an optimized con- 
trol field with Ai = 0.02fttJo, A2 = 0.0035fi,wo, and (j) = 11.44 
rad, while the dots are calculations using the two-mode Flo- 
quet formalism. The dashed curve is a numerical simulation 
with A2 = 0. Here, the system parameters were chosen to 
be aJo/(27r) = 6 GHz and Ns = 4. Other relevant parameters 
are ni/{2Tv) = 86 MHz, Q.2/{2tv) = 15 MHz, uji/{2tv) = 5.785 
GHz, Looi/{2n) = 5.77 GHz, and LJ2/(27r) = LJi2/(27r) = 5.5 
GHz. 



To illustrate this method, we consider a particular ex- 
ample. Fig. [1] shows the result of a numerical simula- 
tion of the time-dependent Schrodinger equation for a 
phase qubit with u!o/{2Tr) = 6 GHz and Ng — 4 subject 
to a control field with Ai = 0.02/iwoi ^2 = 0.0035hujo, 
(j) = 11.44, uji = woi + f^oi/ (2(^01 - W12)), and UJ2 = Wi2- 
The values of A2 and 4> were found by a numerical search 
to optimize the 0^1 transition, providing a significant 
improvement over the A2 — Q dynamics. This search was 
inspired by the general arguments given in Ref. [28| , and 
demonstrates that the use of two frequencies can improve 
the control of this quantum system. 



1 712) with the Floquet Hamiltonian matrix 



{Hq + nifujJi + n2huj2)5n 

+ lA2X{e'^Sn^,n-e, + 6" 



(23) 

where n = {711,77.2}, m = {7711,7^2}, ei = {1,0}, and 
62 = {0,1}. To obtain the state amplitudes, one sums 
over the intermediate photon states 



,{t)= ^ e*("i'^i+"='^2)*(s,77i,772|exp 



"n, 771+62) I Also shown in Fig. [T] is the result of a Floquet cal- 

culation performed by numerically diagonalizing Tip in 
a basis of 22 states, including up to three photons for 
each frequency. Here we provide an analytical approx- 
imation to explain this improved transition. Simula- 
tions suggest that a minimal model for this transition 
involves the states |0,0,0), |1,-1,0), |1,0,-1), |2,-2,0), 
and |2, —1, —1). The Floquet Hamiltonian, in this basis, 
(24) reads 



( ° 


ni/2 


rj2e''^/2 


\ 









l]2e^'^/V2 


r22e-*'^/2 





A 








f7i/V2 





-A2 


V 


ri2e-*'^/\/2 


Q.1/V2 


/ 



(25) 



where Mli — AiXm, MI2 ~ ^22^01, 5 ~ oji — ljqi, A = found through perturbation theory, we find 
ojQi — UJ2, A2 — 2^1 — IV02, and we have let X12 — V^xqi- 
By carefully normalizing and expanding out the terms 

aoit) « cos(f]ii/2) (1 - ^(1 + cos(2</.))^ + 2^e^25t^ 

(26) 
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aiit) 



-isin(f^it/2) ( 1 - 

n 



4A2 n'f 



2 C0S(2(/)) 



and 



02 (t) 



(l + 2e^(^+3^)* - 3e'(^+*)*) ,(27) 



i4, i(A+S)t 



2 f7i 



r^i 2V2A ■ ^ ' 



We see that, in addition to the Rabi oscillation terms 
seen previously, there are terms that oscillate at the fre- 
quencies 5 = uji — woi = ^i/(2A) and A = wqi — '^12- 
The former oscillations are slow, and can typically be ig- 
nored, but the latter oscillations become important near 
the peaks of the Rabi oscillations. One can, in fact, use 
this to optimize the transition. 

At time T = n/Ui, many terms drop out of these am- 
plitudes, and by looking at the leading order terms of a2, 
one finds that it will vanish provided 



"l -in/2 -i(A+35)T 

2A 



(29) 



This condition, in turn, specifies the optimal amplitude 
and the phase of the second microwave drive. Thus, we 
have identified a procedure to optimize the — *■ 1 tran- 
sition by a controlled interference through the Floquet 
state dynamics. Using this value for cj) and fl2, we find 
that the residual error scales as fif /A^, much better than 
the ri^/A^ scaling found for a single frequency transition. 



NUMERICAL OPTIMIZATION 
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FIG. 2: Numerically optimized fl2 as a function of the pri- 
mary Rabi frequency Qi. The dashed curve is the approxi- 
mation 0,2 ~ r2i/(2A) (see text). 
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FIG. 3: Numerically optimized phase as a function of overall 
pulse time T. The dashed curve is the approximation « 

7r/2 + AT (see text). 



The analysis of the preceding section was motivated 
by optimizing numerically the amplitude and phase of 
the second frequency for the — *■ 1 transition. As shown 
above, it was found that by choosing the amplitude and 
phase appropriately, one can obtain significant improve- 
ment in the transition probability using control fields 
with constant amplitude, called square pulses. Here we 
compare the analytical results with the numerically op- 
timized parameters, and show how this approach can be 
used to generate optimized Gaussian pulses [flj . 

First, in Fig. ^ we show the numerically optimized 
^2 — A2Xoi/h as a function of the bare Rabi frequency 
ill — AiXQi/h for a phase qubit with ujQ/(2Tr) — 6 GHz 
and Ns — 4, comparable to recent experiments [l^; other 
parameters can be found in Fig. [TJ For this system, the 
anharmonicity is A/2tt w 260 MHz. We see that the 
analytical result 

f^2,opt « ^ (30) 



provides an excellent approximation for the optimized 
amplitude. Similarly, the optimized phase (f> is plotted 
as a function of the pulse time T in Fig. [3] As with the 
amplitude, the analytical result 

0opt + (31) 

provides an excellent approximation. 

As a further test of this method, we compare the error 
Pe = ^ — pi(T) for this two- frequency pulse with that of 
a single-frequency pulse. This is displayed in Fig. [J] The 
single-frequency pulse is seen to have an error that scales 
as f2^/A^. We see that the two-frequency pulse does a 
significantly better job compared to the single frequency 
pulse, and the error scales as Of/16A*, with oscillations 
of frequency A. 

Finally, using this approach, one can design pulse 
shapes to further optimized the transition. We consider 
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FIG. 4: Error of 1 transition using square pulses. The up- 
per points (squares) are the error of a pulse using a single fre- 
quency with uii = tJoi. The lower points (dots) are the error 
of an optimized two-frequency pulse with uji — loqi + ^1/ (2 A) 
and UJ2 = iLJi2. The upper dashed curve is 3nV(4A^), while 
the lower dashed curve is n*/(16A'') (see text). 
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FIG. 6: Effective spin rotation from — » 2. The probability 
P2{t) to be in state 2 is shown as a function of time. The 
solid curve is a numerical simulation using Ai = O.OlfitJo and 
A2 — AiXoi/xi2 with parameters of Fig. [1] The dashed curve 
is the expected rotation obtained by using the rotating wave 
approximation. 
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FIG. 5: Error of — > 1 transition using Gaussian pulses. 
The upper dotted curve is the error of a pulse using a single 
frequency with u>i = uioi ■ The dashed curve is the error of a 
single-frequency pulse with ioi = o^oi +d2^li/{2A). The lower 
solid curve is the error of an optimized two-frequency pulse 
with ll}2 ~ ^\2 (see text). 



a Gaussian pulse shape 

/(<) = s(i) {Ax cos{LUit) + A2 cos(cj2i + (32) 

with 

sit) = (e-"(i-2t/T)^ _ g-a^ ^ (33) 

where a specifies the shape of thepulse and Na is cho- 
sen such that / s{t)dt = T [ll|, [l^. These pulses are 
optimized using the bare Rabi frequency 

f^i = ^(l + c.(^) (34) 



and drive frequency 

wi=woi + da^^, (35) 

where the dimensionless coefficients and da are var- 
ied to obtain the best transition. These coefficients cor- 
rect for the reduction in Rabi frequency and the ac Stark 
shift discussed previously, and depend on the pulse shape 
parameter a. For Ng — 4 and a — 2, we find that 
Ca=2 = 0.58 and da=2 = 1.245 are required. The error 
using Gaussian pulses with and without the Stark shift 
correction is displayed in Fig. \5\ We see that the single- 
frequency pulse is not effective without these corrections. 
To incorporate the two-frequency pulse, we numerically 
optimize for A2 and (j), and find that it provides a signif- 
icant advantage. Note, however, that the two-frequency 
square pulse outperforms all of the Gaussian pulses for 
small pulse times. 



VI. THREE-STATE OSCILLATIONS 

Recently, multi-frequency control of multiple levels of a 
superconducting circuit has been experimentally demon- 
strated [lo| . This phase qudit was used to emulate spin- 
1 and spin-3/2 quantum systems. Here we look at the 
spin-1 case, and show how the two-mode Floquet theory 
explains the nature of the three-state oscillations at high 
microwave power. 

Figure [S] shows the state 2 probability p2 (t) , when the 
control field is chosen with Ai = O.OlhujQ, ui — wqi, 
L02 = i^\2^ and h£l = AiXm = ^2X12. Using the rotating 
wave approximation, one expects the dynamics to should 
emulate the rotation of a spin-1 system, yielding a prob- 
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ability to be in state 2 of 

While it is expected that there may be Stark shifts, cor- 
rections to the Rabi frequencies, and off-resonant transi- 
tions for this square pulse, a qualitatively new effect is 
seen in the numerical simulation. This is a beating at the 



nF = h 



Note that one 3-level block is isomorphic to a spin oper- 
ator for a spin-f system. 

By performing the lowest order of perturbation theory 
for the coupling between blocks oiTip, one finds that the 
relevant transition amplitude is 

(38) 

This provides an excellent approximation to the beating 
observed in Fig. [51 Note that the perturbation, which is 
proportional to 51/ A, happens to vanish precisely when 
the unperturbed oscillation reaches its maximum {t = 
\/27r/il.) Thus, it is likely that additional effects limit 
this approach to a ^ 2 transition. By extending the 
matrix to 15 states and higher orders in perturbation 
theory, one finds a state 3 population proportional to 



VII. CONCLUSION 

In this paper we have analyzed a set of multi-level ef- 
fects found in superconducting circuits such as the phase 
or transmon qubit when controlled by pulses with two 
microwave frequencies. These involve a combination of 
resonant, off-resonant, and interference effects that are 
of importance for future qubit (or qudit) superconduct- 
ing implementations of quantum information processors. 
Indeed, we have demonstrated that the many-mode Flo- 
quet formalism for multiple frequencies is a useful gener- 
alization of the standard rotating wave approximation. 

First, we used the single- mode formalisim to recover 



frequency A = lui — lu2- 

Using the Floquet formalism, one finds that the 
dominant effect is a coupling between three photon 
blocks of the three-level system, or a total of nine 
states: |0,-1,1), |1,-2,1), |2,-2,0), |0,0,0), |1,-1,0), 
|2,-1,-1>, |0, 1,-1), |1,0,-1), and |2,0,-2). In this 
basis, the effective Hamiltonian is 



(37) 



I 

compact analytical results for corrections to Rabi oscil- 
lations in a three-level system, finding corrections to both 
the resonance condition and the oscillation frequency of 
relative order fi^/A^. These are the ac Stark shift and 
reduction in Rabi frequency seen in existing experiments 
and predicted previously. 

Second, we have shown that simultaneously controlling 
the qubit with two frequencies, one resonant with the — -> 
1 transition (after compensating for the ac Stark shift) 
and the other resonant with the 1^2 transition leads 
to a useful interference effect. This insight was inspired 
by numerical results on square pulses, and found to be in 
excellent agreement. This approach was further extended 
numerically to show that two-frequency Gaussian pulses 
can be developed for the — > 1 transition with significant 
improvements over single- frequency pulses. 

Finally, we have used the Floquet method to explain 
off-resonant couplings that emerge when using the phase 
qudit to emulate a spin system. Here we have found and 
explained a beating that is proportional to r2/A, and 
should be observable in recent experiments, provided it 
is not masked by effects of decoherence. 

In order for these effects to be genuinely useful, one 
would like to extend the optimization of a transition be- 
tween two (or more) states to the optimization of a uni- 
tary operation acting an a superposition of these states. 
Here, however, an interesting difficulty emerges. For the 
single-frequency pulse, with square or Gaussian shapes, 
this is immediate: this control pulse is symmetric under 
time-reversal: /(— i) = f{t)- Consequently, the transi- 
tion from 0^1 and its time reverse from 1 — s- are 
both optimized for a single f(t). For the two- frequency 
pulse, however, f{—t) ^ f{t), and in fact the optimiza- 
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tion developed in Sec. Ill does not perform as well for the 
1^0 transition. Note that this observation sheds some 
light on the two-quadrature approach of Ref. [IB]: the 
class of control pulses advocated there is time-reversal 
symmetric. We expect that combining multiple quadra- 
tures and multiple frequencies will significantly expand 
the control techniques for future experiments. Develop- 
ing simple, accurate, control pulses for multi-level quan- 
tum systems remains a challenging problem for theory 
and experiment. 

APPENDIX 

In this appendix we summarize the perturbative results 
for the cubic oscillator 

i/ = /i^o Qp' + - Ax^y (A.l) 

Here we summarize the Rayleigh-Schrodinger perturba- 
tion expansion for the Hamiltonian H ~ Hq -f- AV^. First, 
one expresses the n-th energy eigenstate, |5'„), in powers 
of A 

OO 

|vI/„) =^A'=|n,fc). (A.2) 

A;=0 

I 



Extending 



This procedure was implemented in Mathematica to 
calculate the eigenvalues up to A^; the A^ expression 
(|A.7[) was found using a more efficient recursion-relation 
method [s^l, and agrees with [34] (provided one lets 
4^ 4^^). These results, when compared with nu- 



In this expansion, \n, 0) is the n-th energy eigenstate of 
Hq, and ]n, k) are the fc-th order perturbative corrections. 
We also expand the energy eigenvalue in powers of A, 



E,, = Y,^''En,k, (A.3) 

fc=0 



where _ffo]'^,0) = Enfi\n,0). Substituting (|A.2p and 
(|A.3p in the eigenvalue equation 



{HQ + XV)\^n) =En\-^„), (A.4) 



equating like powers of A'"', and projecting onto (m, 0] 
allows one to solve for the energies and eigenfunctions: 



Er,,k = {n,0\V\n,k-l) (A.5) 



and 



I 

merical results found by complex scaling jj^, are found 
to be accurate for states n = ^ 3 when A^s > 3. 

In addition to the energy levels, perturbation theory 
also provides expressions for the wavefunctions. For ref- 
erence we list the third-order expression: 



(m, 0]F]n, fc — 1) — En j{m, 0]n, k — j) 
\n.k) = J2 ' "'^^ ^]m,0). (A.6) 

this calculation to A^ one finds 

EnlfujJQ = (n + l/2)-iA2(30n2-|-30n-|-ll)-i^A'^(94n3 + 141n2-|-109n-^31) 
8 32 

--^A^ (115755?!'' -I- 231510^3 -f 278160n2 + 162405n + 39709) 
128 

— ^A^(2282682n5 -1- 5706706^^ -I- 938769071^ -f 837483071^ -|- 4244573n -I- 916705). (A.7) 
2048 



-1-3 -1-6 4-9 

= ]n) + A ^ au(n)\n + fc) + A^ ^ hk(n)\n + k) + ^ Cfc(n)]n -f- k) 



(A.8) 



where \n) = |n, 0) are the eigenstates of the purely harmonic osciUator Hamiltonian, and the nonzero expansion 

coefficients are 



and 
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c_7(n 
c_5(n 
c_3(n 
c_i(n 
c+i(n 

C+3(« 

c+5(n 
c+7(n 
c+9(n 



a_3(n) 
a_i(n) 
a+i(n) 



-^(n(n-l)(n-2))i/^ 



2^2 
3 

2\/2 



,3/2 



(n+ 1)3/2 



6_6(n 

b-2{n 
b+2{n 
6+4 (n 
6+6 (n 

1 



2592\/2 
1 



«+3W = ^((n + l)(n + 2)(n + 3))^/\ 

= ^ (n(n - l)(n - 2)(n - 3)(n - 4)(n - 5))^/^ 

= -^(n(n-l)(n-2)(n-3))^/'(4n-3) 

= ^(n(n-l))^/^(7n2-19n+l) 

= ^((n+l)(n + 2))'/'(7n2 + 33n + 27) 

= ^ ((" + + 2)(n + 3)(n + 4))'/' (4n + 7) 

= lii + + + 3)(n + 4)(n + 5)(n + 6))^/' , 



(n(n - l)(n - 2)(n - 3)(n - 4)(n - 5)(n - 6)(n - 7)(n - 8))^/^ 



192v/2 
1 

960\/2 
1 



1728v^ 
3 



(n(n - l)(n - 2)(n - 3)(n - 4)(n - 5)(n - 6))^/^ (2n - 3) 
(n(n - l)(n - 2)(n - 3)(n - 4))^/^ (SOn^ - 305n + 164) 
(n(n - l)(n - 2))^^^ (48871^ - 217571^ + 4018n - 825) 



64\^ 

3 



ni/2(20n^ + 81n3 + 326^2 + gin + 44) 



64^2 
1 

1728\/2 



(n + l)^/2(20n'' - + 20371^ + 408n + 228) 



((n + l)(n + 2)(n + 3))^^^ (488n3 + 363971^ + 9832n + 7506) 



^ ((n + l)(n + 2)(n + 3)(n + 4)(n + 5))^^^ {80n^ + 465n + 549) 
((n + l)(n + 2)(n + 3)(n + 4)(n + 5)(n + 6)(n + 7))'/' (2n + 5) 



960^2 
1 



192\/2 

1 

2592^2 



((n + l)(n + 2)(n + 3)(n + 4)(n + 5)(n + 6)(n + 7)(n + 8)(n + 9))^^^ , 



(A.9) 
A.IO) 
A.ll) 
A.12) 

A.13) 
A.14) 
A.15) 
A.16) 
A.17) 
A.18) 

A.19) 
A.20) 
A.21) 
A.22) 
A.23) 
A.24) 
A.25) 
A.26) 
A.27) 
A.28) 



One application of these expressions is to calculate the 
(properly normalized) matrix elements of the position op- 



erator 



((*„|M/„)(M/„|vl/„))i/2- 



(A.29) 
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Using the wavefunctions (jA.Sp and matrix elements of 
the previous section, we find 



with corrections of order A^. 



XOfi = 

xo,i = 

Xq,2 ^ 

Xo,3 = 
Xl,l = 
Xl,2 

Xl.3 
X2.2 
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3^ 33 . 

^/2 IIV2 2 
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= 1 



11 
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16 



15, 573 o 
^/6 , 33V6 2 
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21 



A^ 



1113 , 



(A.30) 
(A.31) 
(A.32) 

(A.33) 
(A.34) 
(A.35) 
(A.36) 
(A.37) 

(A.38) 
(A.39) 
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